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We propose the "supersymmetric (SUSY) Yukawa sum rule", a relationship between physical 
masses and mixing angles of the third-generation quarks and squarks. The sum rule follows directly 
from a relation between quark and squark couplings to the Higgs, enforced by SUSY. It is exactly 
this relation that ensures the cancellation of the one-loop quadratic divergence in the Higgs mass 
from the top sector. Testing the sum rule experimentally would thus provide a powerful consistency 
check on SUSY as the solution to the gauge hierarchy problem. While such a test will most likely 
have to await a future next-generation lepton collider, the LHC experiments may be able to make 
significant progress towards this goal. If some of the terms entering the sum rule are measured 
at the LHC, the sum rule can be used (within SUSY framework) to put interesting constraints on 
the other terms, such as the mixing angles among third-generation squarks. We outline how the 
required mass measurements could be performed, and estimate the accuracy that can be achieved 
at the LHC. 



Introduction — Experiments at the Large Hadron Col- 
lider (LHC) have begun probing physics at the TeV 
scale. The primary goal of these experiments is to under- 
stand the mechanism of electroweak symmetry breaking 
(EWSB). The Standard Model (SM) explanation of the 
EWSB, the Higgs mechanism, suffers from the hierarchy 
problem, whose solution very generically requires that 
new physics beyond the SM appear around the TeV scale. 
Among the candidate models of this new physics sug- 
gested by theorists, supersymmetry (SUSY) is perhaps 
the most appealing one. Simple SUSY models are con- 
sistent with experimental data, including gauge coupling 
unification, and fit naturally in string theory. Searches 
for SUSY will be one of the main directions pursued by 
the LHC experiments. 

Assuming that some of the signatures predicted by 
SUSY models are seen at the LHC, the next major task 
for the experiments will be to determine the nature of the 
new particles involved, such as their masses and spins. 
In addition, there is a large number of couplings involv- 
ing the new particles that one can attempt to measure. 
Among those, there is a small set of couplings that is, in 
our opinion, truly special, and deserves special attention. 
These are the couplings that ensure the cancellation of 
the quadratically divergent diagrams contributing to the 
Higgs mass parameter at one loop. Specific relations be- 
tween these couplings and the SM gauge and Yukawa cou- 
plings are required to solve the hierarchy problem, and 
SUSY guarantees that these relations are satisfied. Test- 
ing these relations experimentally would clearly demon- 
strate the role of SUSY in restoring naturalness to the 
EWSB sector. The first goal of this paper is to suggest a 
simple sum rule, which follows unambiguously from one 
such coupling relation, and involves only physically mea- 
surable quantities. The second goal is to outline the set 
of measurements that would need to be performed to test 
this sum rule, and evaluate the prospects for these mea- 
surements at the LHC. While a test of the sum rule will 
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FIG. 1: Quadratically divergent one- loop contribution to the 
Higgs mass parameter in the SM (a), cancelled by scalar su- 
perpartner contributions in a SUSY model (b). 

have to await a next-generation lepton collider, we find 
that the LHC may be able to measure several ingredi- 
ents of the sum rule. Within the framework of SUSY, 
the sum rule can then be used to infer parameters, such 
as stop and sbottom mixing angles, which will be difficult 
or impossible to measure directly. 

SUSY-Yukawa Sum Rule — The strongest coupling of 
the SM Higgs boson is the top Yukawa. At the one- 
loop level, this coupling introduces a quadratically diver- 
gent contribution to the Higgs mass parameter, via the 
diagram in Fig. [l] (a). In SUSY models, this contribu- 
tion is canceled by the diagrams in Fig. [l] (b), with the 
stops, scalar superpartners of the top quark, running in 
the loop. The cancellation relies on the precise relation 
between the top Yukawa and the stop-Higgs quartic cou- 
pling, shown in the figure, which is enforced by SUSY. We 
would like to test this relation experimentally. The most 
direct test, measuring the stop-Higgs quartic vertex, ap- 
pears impossible at the LHC due to extreme smallness of 
all cross sections involving this vertex. However, once the 
Higgs gets a vacuum expectation value (vev), the quartic 
vertex generates a contribution to stop masses, which are 
in principle measurable. The challenge is to isolate this 
term from other contributions to the stop mass matrix. 

SUSY makes two kinds of predictions: (1) it dictates 
a particular particle content (i.e. superpartners to the 
SM fields), and (2) it imposes certain relations between 
couplings of the fields, such as the relation in Fig. [l] 
We want to separate the two, fixing the particle content 
(which we assume could be tested by independent obser- 
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vations), while attempting to test the couphng relation. 

Start with a SUSY-like particle content for the 3^^ gen- 
eration, i.e. a set of scalars with gauge charges 



tL 



(3, 2)i/6, in ^ (3, 1)2/3, ^ (3, l)-i/3- 



(1) 

Leaving the SU{2)l x f/(l)y gauge symmetry unbroken 
and working in the (ti,, ifi)-basis, the only allowed mass 
terms are 
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(in the MSSM these are just the soft masses). Within 
the chosen particle content, we can parameterize EWSB 
mo del- independently by inserting spurious Y^'^. The 
(1,1) entries of the top- and bottom-partner mass ma- 
trices become 
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where v = 246 GeV. Let us define an observable 
1 
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where the top-partner eigenmasses rriti < mt2, the 
bottom-partner eigenmasses 77151 < 77152, and the mixing 
angles Ot and ^5 are all, in principle, measurable. (We 
use the notation Ct^b = cos^t^5, ^^^5 = sin^^^5.) Writing 
the top-partner mass matrix in terms of these quantities: 

(similarly for M-') and canceling the soft mass M£ by 



evaluating (M|)ii — (M|)ii, we obtain 
T = FA - . 
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In other words, T probes the spurious only. Note, how- 
ever, that Eq. ([6| will receive non-trivial corrections be- 
yond the tree level, since T is defined in terms of physical 
(pole) masses, while in the above derivation all masses are 
evaluated at the same scale. 

At tree level, SUSY makes a definite prediction for T. 
Using the standard sfermion tree-level mass matrices (see 
e.g. [1 ) and neglecting fiavor mixing, we obtain 
1 
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Here the hats denote tree- level (or "bare") masses. The 
numerical values assume the renormalization scale Q = 
600 GeV (so that i.e. rht ~ 153 GeV), but do not depend 
strongly on the precise value of Q. This prediction, which 
we call the SUSY- Yukawa sum rule^ relies on the same 
relation between the fermion and scalar Higgs couplings 
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FIG. 2: Distribution of T for a SuSpect random scan of 
pMSSM parameter space. Scanning range was tan /3 G (5, 40); 
Ma, Ml e (100,500) GeV; M2, M3,\fi\, Mql, MtR, MbR G 
(Ml + 50 GeV, 2 TeV); |At|, < 1.5 TeV; random 

sign(/i). EWSB, neutralino LSP, and experimental con- 
straints (?7T.i/, Ap, b 57, a^, m~± bounds) were enforced. 



which leads to the cancelation in Fig. [T] Measuring T 
would therefore provide a powerful, if somewhat indirect 
method of testing whether it is SUSY that solves the hi- 
erarchy problem. (This argument is conceptually similar 
to the tests of the Little Higgs cancellation mechanism, 
proposed in [2j. Earlier examples of SUSY sum rules, 
devised within the mSUGRA framework, can be found 
in i.) 

Radiative corrections to the SUSY prediction for T 
can be important, since the sum rule typically involves 
a rather delicate cancellation between stop and sbottom 
mass terms. The full analytical expressions for the radia- 
tive corrections to superpartner masses within the MSSM 
can be found in [1 , and a convenient numerical imple- 
mentation is provided by the SuSpect package [4]. The 
corrections depend on a large number of MSSM parame- 
ters. To estimate their effect on T, we conducted several 
scans of the MSSM parameter space using SuSpect. We 
did not assume a particular model of SUSY breaking, but 
allowed the weak-scale soft terms to vary independently. 
A representative result for the distribution of T is shown 
in Fig. [2] (As usual, the reader must exercise caution 
in interpreting this plot, since it necessarily reflects our 
sampling bias of parameter space.) It shows that radia- 
tive corrections can change the value of T significantly 
from its tree level prediction However, a measure- 
ment of |T| > 0(1) would disfavor TeV-scale SUSY as 
the solution to the hierarchy problem. It should be 
noted that in a generic theory with the particle con- 
tent of Eq. ([1]), the scalar-Higgs quartic couplings are 
only constrained by perturbativity, leading to the possi- 
ble range of — IGtt^ < T < IGtt^. Moreover, if some of 
the parameters in the sum rule are misidentified, an even 
broader range is possible. For example, if the mixing an- 
gle measurements were off by 7r/2, the right-hand side 
of Eq. Q would contain the right-right elements of the 
squark mass matrices, which are of course independent 



3 



for stop and sbottom, so any value of T is in principle 
possible. Thus, even with radiative corrections included, 
the SUSY- Yukawa sum rule presents a useful and non- 
trivial consistency check on SUSY. 

It is also interesting to ask if the sum rule can be 
used as a tool for model discrimination. Recently, several 
SUSY "look-alikes", i.e. models whose LHC signatures 
are similar to SUSY but arise from completely different 
underlying physics, have been studied. The most stud- 
ied examples are universal extra dimensions (UED) |5] 
and little Higgs with T-parity (LHT) [6^ models. These 
models contain particles with the quantum numbers of 
Eq. ([T]), but instead of scalars, they are spin- 1/2 fermions. 
(The minimal LHT model does not contain a hR coun- 
terpart; however, such a particle can easily be added.) 
This leads to a different Higgs coupling structure: for ex- 
ample, the 4-point coupling in Fig. 1 (b) does not exist, 
at renormalizable level, in these theories. As a result, 
UED and LHT predictions for T are generically different 
from SUSY, at least at the tree level. As an example, the 
tree-level prediction of the minimal LHT model is 
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where m^^ and ttiah the masses of the heavy, T- 
odd partners of the left-handed b quarks and the hyper- 
charge gauge boson, respectively. In contrast to SUSY, 
T is always negative at tree level in the LHT; for typi- 
cal parameter values T ^ —0.5. Unfortunately, radiative 
corrections can shift T in SUSY significantly, including 
changing the sign, as can be seen in Fig. |2] Presum- 
ably, the LHT prediction will also receive important loop 
corrections, although they have not yet been calculated. 
Depending on the resulting ranges and on the measured 
value of T, the measurement may be interpreted as sup- 
porting one or the other model, but it seems unlikely 
that a sharp model-discriminating statement could be 
made. On the other hand, one should keep in mind that a 
measurement of parameters not directly entering the sum 
rule (such as the gluino mass) would generally shrink the 
range of possible T values in each model by constraining 
the possible radiative corrections, improving the model- 
discriminating power of this observable. 

Measuring all the ingredients of T is very difficult at 
a hadron collider, and the determination of the complete 
3^^-generation sfermion spectrum and mixing angles will 
most likely have to be performed at a future lepton ma- 
chine. However, for favorable MSSM parameters, some 
progress can be made at the LHC. In particular, if some 
of the ingredients of the sum rule can be measured, and 
the sum rule is assumed to be valid, it can be used to 
put interesting constraints on the remaining ingredients. 
The easiest terms to measure at the LHC are the masses 
of the lightest stop and sbottom squarks. To understand 
the implications of such a measurement, let us rewrite T 



H2 



Hi 



s2 



^62 



(9) 

Assuming that the SUSY framework is correct, a mea- 
surement of together with the sum rule can be used 
to constrain the third-generation mixing angles, even if 
nothing is known about the masses of the heavier su- 
perpartners ?2 and 62- This is illustrated by the scatter 
plots in Fig. [sj If is small, then either both ii and 
bi must be mostly left-handed so that ATt, 5 is small, 
or the two AT's must precisely cancel each other. (Ob- 
viously, the second possibility is less likely, as reflected 
in the distribution of points in Fig. [3] (b).) A large and 
negative would require a right-handed ti, whereas a 
large and positive requires a right-handed bi. Thus, 
mass measurements together with the sum rule can pro- 
vide non-trivial information on the mixing angles, which 
are difficult or impossible to measure directly at the LHC. 
(For some proposals for measuring the stop mixing angle, 
see Refs. [II [H].) 

Prospects at the LHC: a Case Study — The MSSM pa- 
rameter point we will consider is defined by the following 
weak-scale inputs (from here on all masses in GeV unless 
otherwise noted): 
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with all other A-terms zero and all other sfermion soft 
masses set at 1 TeV. The relevant spectrum (calculated 
with SuSpect) is 
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At this benchmark point, T = 0.423, and = 0.350. 
We will show below that the LHC can measure rather 
accurately. 

To measure the ti and bi masses, we propose to use 
kinematic edges, the classical Mt2 variable [9 , and re- 
cently proposed "subsystem- Mt2" variables \iO^ to ana- 
lyze the two processes 

(I) g ^ bib ^ bbxi via gluino pair production, 

(II) ii txi via stop pair production 

(where we omit antiparticle indices). For our bench- 
mark point each of the above decays has 100% branch- 
ing fraction, completely eliminating irreducible SUSY 
backgrounds to the measurements discussed below. The 
process (I) yields the ^,61, and Xi masses, and the pro- 
cess (II) provides rriti. Below, we briefly outline these 
measurements, and estimate their accuracy; the details 
of this analysis will be presented in [11 . 

We ignore issues related to hadronization and ISR by 
performing the analysis at leading order in as and at 
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FIG. 3: Scatter plot of pMSSM parameter points produced by the SuSpect scan from Fig.|2] showing the correlations between 
the stop and sbottom mixing angles for different ranges of T^ Each 0.005 x 0.005 bin is colored according to the number of scan 
points contained in it, with hot (bright) and cold (dark) colors indicating high and low scan point density, and unpopulated bins 
left uncolored. These correlations are a direct consequence of the SUSY- Yukawa Sum Rule, and any measurement of ^ 
provides valuable information about the sbottom mixing angle. 



parton level. We use MadGraph/MadEvent ( MGME) 
package [12 to simulate gluino and stop production, 
and BRIDGE [13] to simulate decays. We use the 
CTEQ611 [14] parton distribution functions through- 
out, with the MGME default (pT-dependent) factor- 
ization/renormalization scale choice. To roughly model 
detector response to jets and electrons, we introduce a 
Gaussian smearing of their energies according to [15] 



50% 



= ©3% 



10% 



GeV 



GeV 



:©0.7%. (10) 



(I) Measuring the bi^g^Xi 'rnasses — We study gluino 
pair production with subsequent decay into 4b + 2xi 
at the LHC with = 14 TeV and 10 fb"^ of inte- 
grated luminosity. The selection cuts are as follows: (a) 
pr > 200 GeV, (b) exactly 4 tagged b-jets, (c) p?)^^ > 
100 GeV, (d) > 40 GeV, (e) \r]\ < 2.5, AR > 0.4. 

The gluino pair production cross section is agg ^11.6 pb. 
We assumed a 6-tag efficiency of 0.6 and b-mistag rates 
for C-, r-, and light quark/gluon jets of 0.1, 0.1 and 0.01, 
respectively, leaving about 1.5 pb of fully 6-tagged sig- 
nal. The other kinematic cuts (a, c-e) have an efficiency 
of 32%, yielding 480 fb, or about 4800 signal events at 
10 fb"^ 

We computed the cross sections of the two main SM 
background processes, 4j + Z with Z ^ vu, and tt with 
one or both tops decaying leptonically. The cross sec- 
tions, including efficiencies of the cuts (a-e), are < 10 fb 
and 25 fb, respectively. Thus, we conclude that the SM 
backgrounds can be effectively eliminated by cuts, and 
do not take them into account further in the mass deter- 
mination analysis. 

The main background for mass determination comes 
from combinatorics. Consider the dijet invariant mass 



M55. If both 6's come from the same decay chain, the 
distribution has a kinematic edge at 
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For each event, there are three possible ways to assign 
4 6's to two decay chains, and the M55 distributions of 
the wrong combinations extend well beyond M^^^. If all 
combinations are included, the edge is washed out. We 
find that the combinatoric background can be reduced 
with simple cuts: very generally, the directions of jets 
from the same decay chain should be correlated, and the 
pairings with the largest invariant masses are likely to be 
incorrect. Denoting the two 6's assigned to each decay 
chain as (1,2) and (3,4) respectively, we drop the com- 
bination with the largest Max[Mi2,M34] in each event, 
and require Max[Ai?i2, Ai?34] < 2.5. The resulting dis- 
tribution shows a clear edge. We fit to it with a simple 
trial-PDF, the linear kink function, which we will use 
throughout this analysis: 



Xi K X2 

An unbinned maximum-likelihood fit reliably finds the 
edge position K, yielding a measurement of the kine- 
matic edge position M^b^^H^ = (395 ± 5) GeV. This is 
quite close to the correct value, Eq. ( [TT] ), but the use of 
the simple linear fit function clearly does introduce a sys- 
tematic error into the edge measurement. To account for 
this effect, we will simply assume a systematic error of 
3 times the statistical error for each edge measurement; 
this is sufficient to bring across the main points of our 
analysis. More sophisticated methods for kinematic edge 
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Mt2^'°(0) without combinatorial error reduction 
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Mt2 (0) with known decay chain assignment 
K = (310.7 ± 3.7) GeV 
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FIG. 4: M^\^{0) distributions. The analytical prediction for the edge position is 320.9 GeV. We emphasize that even though 
we show the linear kink fits only over a certain range, K depends very little on the fit domain. 



extraction exist in the literature (e.g. [7 ), and would be 
used in practice. 

The position of the kinematic edge provides one func- 
tion of the three unknown masses; two more are required 
to solve for the spectrum. These can be obtained from the 
endpoints of distributions of events in MT2-subsystem 
variables [10^ M^f{0) and M|^°(0), predicted to be at 
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TABLE I: Mass measurements (all in GeV), assuming Gaus- 
sian edge measurement uncertainties. We imposed the lower 
bound m^o > 45 GeV, which generically follows from the 
LEP invisible Z decay width measurement fT7\ . 



M|f (0)"'''" = m~g- mlo/m~g = 506.7 GeV. 



Of the several possible Mt2 variables for this system, 
these two show the clearest edges, allowing precise mass 
determination; the complete analysis of all Mt2 variables 
will be presented in pT] . 

To calculate M^\^ for each event, we must divide the 
four 6's into an upstream and a downstream pair, giving 
6 possible combinations. Fig. [4] (a) shows the complete 
M^^{0) distribution; the edge is completely washed out. 
It turns out that of the 5 possible wrong pairings, the two 
where 6's from the same decay chain are put into up- and 
down-stream pairs are the most problematic, since their 
distributions extend significantly beyond the edge. 
Based on this observation, we developed two techniques 
to reduce the combinatorial error. Firstly, for each event 
we can simply drop the two largest Mt^s. The corre- 
sponding distribution is shown in Fig. [4] (b) . Secondly, 
we can use our measurement of the kinematic edge. For 
each event there are three possible ways to assign the 4 
6's to two decay chains. For some events (about 30% in 
our sample) we find that for two of these combinations, at 
least one same-chain invariant mass is larger than M^^^, 
whereas for the other combination both same-chain in- 
variant masses are smaller - this combination must be 
the correct one. Using only those events and keeping 
only the correct decay chain assignments, we obtain the 
distribution of M|,2^(0) shown in Fig. [i] (c). We per- 
formed linear kink fits on the distributions in Fig. [4] (b) 
and (c), and found that they are in agreement, indicat- 
ing the robustness of our approach. Combining the two 
fits yields M^\^{0)2Zs = (314.0 ± 4.6) GeV. We used a 



(12) similar method to extract the edge, and obtained 

^If (O)meas = (^^^.l ± 4.8) GcV. As for the kinematic 
edge, the linear fit function works rather well, but it does 
introduce some systematic error into the edge measure- 
ments, which we again model by inflating the error bars 
by a factor of 3. To summarize, the measured edges are: 



max 
^meas 



M|f(0) 



meas 
max 



(395 ± 15) GeV , 
(314 ± 14) GeV , 
(492 ± 14) GeV . 



(13) 



Each of these edges defines a subvolume of 
(m^, m^o, m^J-space, which yields the mass mea- 
surements given in Table [ij 

(II) Measuring the ii-mass — We simulate pp 
iiil ti^ 2xi for 100 fb"^ integrated luminosity. The 
signal production cross section is 2 pb. The dominant 
irreducible background is {Z vv)tt with (Jbg = 135 
fb. Following [16 , we demand two fully reconstructed 
hadronic tops in each event, in order to use the classical 
Mt2 variable [9]. Our signal cuts are (a) exactly 2 tagged 
b-jets and at least 4 other jets with pt > 30 GeV and 
|7^| > 2.5 (b) lepton veto (c) AR > 0.4 between all the b- 
and light jets (d) pr > 100 GeV (e) Ht > 500 GeV (e) 
^max ^ j^QQ require 4j to reconstruct to two l^'s 

with a mass window of (60, 100) GeV and the two VK's to 
reconstruct with the two 6's to two tops with a mass win- 
dow of (140,200) GeV. After cuts we are left with 1481 
signal and 105 background events. Plotting the classical 
Mt2 distribution we see a clear edge, and using the linear 
kink fit trial PDF with error scaling yields 

Mt2(0)SS, = (340±4) GeV. (14) 
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Compare this to the analytical prediction [T8| 
Mt2(0)""^^ = 336.7 GeV. Combined with the m^o 
measurement from (I), this yields the stop mass m^i, see 
Table [T[ Taking into account all correlations, we find: 

'^'mee. = ^ Kl " = 0.525t°;?° , (15) 

in good agreement with the theoretical value = 0.350. 
As explained above, a measurement of does not by 
itself provide a consistency check of SUSY, or help in 
discriminating it from other models. However, if the 
SUSY- Yukawa sum rule is assumed to be valid, this mea- 
surement can be used to place a constraint on the 3rd 
generation squark mixing. The measurement in Eq. ([T5| 
corresponds to the range of assumed in Fig. [s] (b). 
Thus, even without using information from any other 
measurements, one could conclude that, most likely, the 
stop and sbottom mixing angles are rather small and the 
observed light stop and sbottom states are mostly left- 
handed (although right-handed light states, with an ac- 
cidental cancellation of AT5 and AT^, would remain as 
a logical possibility at this point). 

Discussion and Conclusions — In this paper we pro- 
posed the SUSY Yukawa sum rule with direct connection 
to the cancelation of quadratic Higgs mass divergence, 
and introduce an observable T that can be used to test it. 
This constitutes a significant check on TeV-scale SUSY 
as the solution to the hierarchy problem. While full mea- 
surement of T will have to be left to a future lepton ma- 
chine, we have demonstrated that progress could already 
be made at the LHC. In particular, we showed that, for 
the MSSM benchmark point we chose, two masses enter- 
ing the sum rule, m^i, can be measured. Given these 
measurements, one could then use the sum rule (within 
the SUSY framework) to put interesting constraints on 
other parameters, such as third- generation squark mix- 
ing angles, whose direct measurement would be difficult 
or impossible. 

In the course of the analysis we developed new tech- 
niques for reducing combinatorial background for Mt2- 
measurements, allowing for complete mass determination 
of and Xi- At this point, we performed the 

analysis at the parton level, with only a crude Gaus- 
sian smearing to account for detector effects. It is im- 
portant to confirm the proposed techniques with more 
detailed simulations including initial and final state radi- 
ation, showering and fragmentation, and better detector 
modeling. Results of a study including some of these ef- 
fects will be presented in Ref. [11 . In the future, it will 
also be interesting to assess the abilities of the LHC to 
test the sum rule (fully or partially) in the MSSM param- 
eter regions with spectra different from our benchmark 
point, as well as to study in detail how the sum rule tests 
can be completed at a future lepton collider. 



Acknowledgments — We are grateful to James Alexan- 
der and Konstantin Matchev for useful discussions. This 
work is supported by the U.S. National Science Founda- 
tion through grant PHY-0757868 and CAREER award 
PHY-0844667. MB thanks the Galileo Galilei Institute 
for Theoretical Physics for the hospitality and the INFN 
for partial support during the completion of this work. 



[1] D. M. Pierce et al. Nucl. Phys. B 491, 3 (1997) 
arXiv:hep-ph/9606211] . 

[2] G. Burdman, M. Perelstein and A. Pierce, Phys . 
Rev. Lett. 90, 241802 (2003) arXiv:hep-ph/0212228; ; 
M. Perelstein, M. E. Peskin and A. Pierce, Phys. Rev. 
D 69, 075002 (2004) arXiv:hep-ph/0310039 . 

[3] S. P M artin a nd P Ramond, Phys. Rev. D 48, 5365 
(1993) |arXiv:he p-ph/9306314 . 

[4] A. Djouadi, J. L. Kneur and G. Moultaka, Comput. Phys. 
Commun. 176, 426 (2007) arXiv:hep-ph/0211331 . 

[5] T. Appelquist, H. C. Cheng and B. A. Dobrescu, 
Phys. Rev. D 64, 035002 (2001) arXiv:hep-ph/0012100 ; 
H. C. Cheng, K. T. Matchev and M. Schmaltz, Phys. 
Rev. D 66, 056006 (2002) arXiv:hep-ph/02053l4 . 

[6] N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nel- 
son, JHEP 0207, 034 (2002) arXiv:hep-ph/0206021 ; 
H. C. Cheng and 1. Low, JHEP 0309, 051 (2003) 
[aFXiv:hep-ph/0308199 . For a review and further refer- 
ences, see M. Perelstein, Prog. Part. Nucl. Phys. 58, 247 
(2007) arXiv:he p-ph/0512128 . 

[7] J. Hisano, K. Kawagoe and M. M. Nojiri, Phys. Rev. D 
68, 035007 (2003) arXiv:hep-ph/03042l4 . 

[8] M. Perelstein and A. Weiler, JHEP 0903, 141 (2009) 
arXiv:0811.1024 [hep-ph]]; K. Rolbiecki, J. Tattersall 
and G. Moortgat-Pick, arXiv:0909.3196 [hep-ph]. 

[9] A. Barr , C. Les ter and P. Stephens, J. Phys. G 29, 2343 
(2003) |arXiv:h ep-ph/0304226 . 
[10] M. Burns, K. Kong, K. T. Matchev and M. Park, JHEP 

0903, 143 (2009) arXiv:0810.5576 [hep-ph]]. 
[11] M. Blanke, D. Curtin, M. Perelstein, in preparation. 
[12] J. Alwall et a/., JHEP 0709, 028 (2007) arXiv: 0706.2334j 

[hep-ph] ] . 

[13] P. Meade and M. Reece, arXiv:hep-ph/0703031> 

[14] J. Pumplin et. al., JHEP 0207, 012 (2002) ar Xivihep- 
ph/0201195 . 

[15] ATLAS Collaboration TDR Vol. 1. CERN-LHCC-99-14, 
ATLAS Collaboration TDR Vol. 2. CERN-LHCC-99- 
15; G. L. Bayatian et al [CMS Collaboration], CERN- 
LHCC-2006-001; G. L. Bayatian et al [CMS Collabora- 
tion], J. Phys. G 34, 995 (2007). 

[16] P. Meade and M. Reece, Phys. Rev. D 74, 015010 (2006) 
arXiv:hep-ph/0601124 . 

[17] M. Acciarri et al [L3 Collaboration], Phys. Lett. B 431, 
199 (1998); R. Akers et al [OPAL Collaboration], Z. 
Phys. C 65, 47 (1995); D. Buskulic et al [ALEPH Col- 
laboration], Phys. Lett. B 313, 520 (1993). 

[18] W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, JHEP 
0802, 035 (2008) |arXiv:071 1.45 26 [hep-ph]]. 



